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Abstract.  For  finite  deformations  of  nonlinear  viscoelastic  solids,  the  speed  of  propagation  of  acceleration 
waves  (i.e.,  ramp  waves)  generally  depends  not  only  on  the  current  state  of  strain  at  the  wave  front  but  also  on 
the  prior  strain  history.  Consequently,  explicit  formulas  for  the  wave  speed  can  be  quite  complicated.  Simple 
formulas  for  the  wave  speed  do  exist  for  special  classes  of  materials  and/or  special  deformation  histories,  and 
in  this  regard  we  consider  one-dimensional  motions  of  viscoelastic  solids  governed  by  single  integral  laws. 
Some  of  the  relations  obtained  are  universal  in  the  sense  that  they  hold  for  all  materials  in  a  given  class  and 
do  not  explicitly  involve  the  relaxation  kernel  function  in  the  hereditary  integral  defining  these  materials. 


INTRODUCTION 

We  consider  the  speed  of  propagation  of  accelera¬ 
tion  waves  in  viscoelastic  solids  undergoing  uniaxial 
strain,  as  would  occur  in  a  normal  plate  impact  exper¬ 
iment  prior  to  the  arrival  of  lateral  release  waves.  The 
front  of  a  ramp  wave  is  an  example  of  an  acceleration 
wave.1  Expansive  (unloading)  ramp  waves  can  be 
generated  by  reflection  of  a  shock  wave  from  a  free 
surface  or  a  lower  impedance  material  [1].  Compres¬ 
sive  ramp  waves  can  be  generated  by  use  of  fused 
silica  buffer  plates  or  by  graded  density  impactors 
[1,  3]  and  also  by  fast  pulsed  power  techniques  [3], 

For  nonlinear  viscoelastic  solids,  the  acceleration 
wave  speed  U  is  generally  a  nonlinear  function  of 
the  current  strain  as  well  as  the  past  strain  history 
at  the  material  point  instantaneously  situated  on  the 
wave  front.  In  particular,  for  viscoelastic  materials 
governed  by  single  or  multiple  integral  laws,  the  de¬ 
pendence  of  the  wave  speed  on  the  strain  history  typ¬ 
ically  involves  an  explicit  dependence  on  the  relax¬ 
ation  kernel  function(s)  in  the  hereditary  integral(s).2 

However,  for  special  classes  of  viscoelastic  solids 
and/or  special  strain  histories,  simple  explicit  formu¬ 


las  for  the  wave  speed  exist  in  terms  of  quantities 
which  have  a  direct  physical  interpretation.  An  ex¬ 
ample  of  such  a  relation  was  given  by  Nunziato  et 
al.  [1]  for  an  acceleration  wave  propagating  into  a 
deformed  region  in  equilibrium  in  a  finite  linear  vis¬ 
coelastic  solid.3  They  showed  that 


nrr2_d0E  ,  l-6t  °l  (£1 )  —  °E  (£1 ) 

P°  de,  i  -  ie,  e 


where  cr  and  are  the  instantaneous  and  equilib¬ 
rium  elastic  response  functions,  and  £j  is  the  equilib¬ 
rium  uniaxial  strain  ahead  of  the  wave.4  This  relation 
is  universal  in  the  sense  that  it  holds  for  all  materials 
in  the  indicated  class  and  does  not  explicitly  involve 
the  relaxation  kernel  function  G'  in  the  hereditary  in¬ 
tegral  (17)  defining  the  materials  in  this  class. 

The  paper  begins  with  a  discussion  of  the  gen¬ 
eral,  one-dimensional,  nonlinear,  single  integral  law 
for  viscoelastic  response.  This  is  followed  by  some 
general  results  on  acceleration  wave  speeds  in  such 
materials.  These  results  are  used  to  derive  relations 
analogous  to  (1)  but  for  more  general  classes  of  non¬ 
linear  viscoelastic  solids  and  non-equilibrium  condi¬ 
tions  ahead  of  the  wave. 


1  More  precisely,  an  acceleration  wave  is  a  propagating  singular 
surface  across  which  the  stress,  strain  and  particle  velocity  are 
continuous  but  their  spatial  gradients  and  time  derivatives  suffer 
jump  discontinuities  [1,  2], 

2  See  eq.  (13)  below  for  the  general  single  integral  case. 


3  This  class  of  nonlinear  viscoelastic  materials  is  defined  by  equa¬ 
tions  (17)  and  (20)  below. 

4  Precise  definitions  of  all  terms  are  given  below. 


CP706,  Shock  Compression  of  Condensed  Matter  -  2003,  edited  by  M.  D.  Furnish,  Y.  M.  Gupta,  and  J.  W.  Forbes 
2004  American  Institute  of  Physics  0-7354-0181-0/04/$22.00 
263 


SINGLE  INTEGRAL  LAWS 


Let  F  denote  the  deformation  gradient  relative 
to  a  fixed,  unstressed  reference  state,  let  J  -  detF, 
and  let  T  denote  the  Cauchy  stress  tensor.  Then 
E  =  7T(F~1)r  and  S  =  7F-1T(F-1)r  are  the  1st 
and  2nd  Piola-Kirchhoff  stress  tensors.  Introduce  a 
Cartesian  coordinate  system  with  the  1-axis  a  sym¬ 
metry  axis  of  the  material,  and  consider  a  time- 
dependent  uniaxial  strain  along  this  axis.  The  normal 
component  of  stress,  taken  positive  in  compression, 
is 

=  —  =  —  Eji  =  —  Aj  ,  (2) 

where  Xl  =  Fn  is  the  (principal)  stretch  along  the  1- 
axis.  Let  e  denote  the  nominal  measure  of  uniaxial 
strain  taken  positive  in  compression:  £  =  1  —  Xl . 

We  consider  viscoelastic  materials  governed  by  a 
nonlinear  single  integral  law,  the  most  general  one¬ 
dimensional  form  of  which  is 

poo 

a(t)  =  oi(e(t))+  J  i f'(e(t),e(t-s),s)ds.  (3) 

Here  (S'  denotes  the  partial  derivative  of  f!  with 
respect  to  its  third  (or  temporal)  argument  s : 

&'(eve2,s)  =  ^  &(eve2,s).  (4) 

Fading  memory  requires  that  the  relaxation  kernel 
Sf ’ (£j ,  £2 ,  s)  decay  to  zero  sufficiently  rapidly  as  s  — > 
°°.  There  is  some  non-uniqueness  in  the  functions  a 
and  .  This  may  be  removed  by  the  assumptions 

^'(£j,0,s)=0  and  &(e1,e2,0)  =  <7;(£j) .  (5) 


Indeed,  since  rather  than  itself  appears  in  the 
constitutive  relation  (3),  we  are  free  to  choose  the 
initial  value  ^(£j,£,,0).  The  choice  (5)0  simplifies 
the  results  below.  By  (5)j  the  upper  limit  °°  in  (3) 
may  be  replaced  with  t  whenever  e(t)  =0  for  T  <  0. 
Condition  (5)1  implies  that  <7  is  the  instantaneous 
elastic  response  function,  i.e.,  cr(f)  =  <r  (£j)  for  the 
jump  strain  history 


£(T) 


£j  ,  if  T  =  t ; 
0 ,  if  T  <  t ; 


(6) 


then  <7  (0)  =  0.  Also  note  that  (5)  implies 

Sf(e1,0,s)=Sf(e1,0)0)  =  a/(e1).  (7) 


Observe  that  (6)  is  the  strain  history  experienced 
by  a  point  which,  at  the  instant  t,  lies  on  the  front  of 
a  shock  wave  propagating  into  an  unstrained  region. 
Thus  the  instantaneous  elastic  response  function  may 
be  inferred  from  measurements  of  the  stress  jump 
across  shocks  in  undeformed  materials  [1], 

For  any  material  with  fading  memory,  let  cr^  (ex ,  t) 
denote  the  stress  at  time  t  >  0  for  the  stress  relaxation 
test 


£(T) 


£j  ,  if  T  >  0 ; 
0 ,  if  T  <  0 . 


(8) 


Then  cr  is  called  the  stress  relaxation  function  (cf. 
[4]),  and 


o-r(£i,0)  =  c7;(£i),  <r<f(e1,oo)  =  CT£(e1),  (9) 

where  a  is  the  equilibrium  elastic  response  func¬ 
tion.  From  (3)— (5)  it  follows  that  for  t  >  0, 

y{evEvf)  =  oR(evt),  (10) 

and  that  (£t ,  e2,t)  is  the  stress  at  time  t  >  0  for  the 
strain  history 

{£,  ,  if  T  =  / ; 

£2,  if  0  <  T  <  f ;  (11) 

0,  if  T  <  0 . 

The  1-D  linear  theory  of  viscoelasticity,  namely 

poo 

a(t)  =  G(0)e{t)+  G’{s)e{t-s)ds,  (12) 

Jo 


is  the  special  case  of  (3)  with  Sf(£j,£,,s)  =  G(s)e2. 
Here  G'(s)  =  dG/ds ;  G(0)  and  G(°°)  are  the  instan¬ 
taneous  and  equilibrium  elastic  moduli;  cr  (e(f))  = 
G(0)e(t )  is  the  instantaneous  elastic  response;  and 
<7R(£1;f)  =  G(t)e j.  In  this  case  G  rather  than  <jr  is 
usually  referred  to  as  the  stress  relaxation  function. 


ACCELERATION  WAVE  SPEED 

Let  po  and  p  be  the  densities  in  the  undeformed 
and  deformed  states.  Let  U(t)  denote  the  referential 
or  Lagrangean  acceleration  wave  speed  (measured 
with  respect  to  distance  in  the  undeformed  reference 
state).  For  a  general  class  of  materials  with  fading 
memory,  Coleman  et  al.  [2]  showed  that  p.fj2{t)  is 
given  by  the  derivative  of  the  stress  response  func¬ 
tional  with  respect  to  the  current  strain  e(t),  holding 
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the  past  strain  history  fixed.  When  applied  to  the  sin¬ 
gle  integral  law  (3),  this  yields  the  formula 

P0U2(t)  =  a'(e{t)) 

r°°  (13) 

+  J  <91^',(e(f),e(f  —  s),s)ds, 

where  <7  '(e)  =  <7  (e)  and  d {S'  denotes  the  partial 

derivative  of  CS'  with  respect  to  its  first  argument. 
The  Eulerian  wave  speed  (measured  with  respect  to 
distance  in  the  deformed  state)  is  given  by  =  XlU . 
And  since  po  =  Ajp,  we  also  have  pfy2  =  Aj  •  pnU2. 

For  the  linear  theory  (12),  we  see  that  (13)  re¬ 
duces  to  poU2  =  G( 0),  and  hence  we  recover  the 
well-known  result  that  the  wave  speed  is  a  constant, 
independent  of  the  deformation  ahead  of  the  wave. 
However,  for  the  nonlinear  theory,  (13)  implies  that 
the  acceleration  wave  speed  U  ( t )  is  generally  a  com¬ 
plicated  function  not  only  of  the  current  strain  e(t) 
but  also  of  the  past  strain  history  at  the  material  point 
instantaneously  situated  on  the  wavefront. 

If  the  material  ahead  of  the  wave  was  initially 
undeformed  and  subjected  to  a  step  in  uniaxial  strain 
of  amount  £j  at  time  zero  (i.e.,  the  stress  relaxation 
test  (8)),  then  (13)  simplifies  to 

P0U2(t)  =  dl&(evel,t) 

=  d1oR(evt)-d2&(evel,t) 

for?  >  0,  where  (14),  follows  from  (10). 

More  generally,  for  the  strain  history  (11)  ahead  of 
the  wave,  the  acceleration  wave  speed  at  the  instant  t 
is  given  by 

P0U2(t)  =  d^(eve2,t).  (15) 

Actually,  this  statement  requires  some  qualification 
since  £  is  not  continuous  at  the  instant  t.  Let  t  >  0  be 
fixed,  and  consider  the  strain  history 


FINITE  LINEAR  VISCOELASTICITY 

Now  we  consider  single  integral  laws  of  the  form 

a(f)  =  cr/(e(0) 

r°  (17) 

+  Jo  G'(e^,s^  '/(e(7-s))c/'s’ 

where  G'  (£1’S)  “  ds  G(£j,s).  This  is  the  special  case 
of  (3)  with 

^'(£1,£2,i)=G,(£1,s)  •/(£,).  (18) 

To  satisfy  (5)2  we  require  that 

&(£i , S2,s)  =  /(£,)  [G(£j , s)  -  G(£j , 0)] 

+  (ei )  ■ 

Here  /  is  interpreted  as  a  strain  measure,  so  that 
/( 0)  =  0  and  /'( 0)  =  1 ,  in  which  case  (5)j  is  sat¬ 
isfied.  Nunziato  et  al.  [1]  considered  the  special  case 
of  (17)  with  /  given  by 

f(e)  =  e-\e2.  (20) 

This  case  arises  from  the  3-D  single  integral  law 

poo 

S(0  =  ^(E(0)+  /  G,(E(f),«)[E(f-j)]rfi, 

(21) 

where  E  =  \  (FTF  —  I)  is  the  Green  strain  tensor 
and  G' (E is  a  fourth  order  tensor.  This  class 
of  viscoelastic  materials  was  termed  finite  linear  by 
Coleman  and  Noll  [5].  The  quadratic  term  in  (20)  re¬ 
sults  from  conversion  to  the  nominal  strain  measure 
£.  Replacing  E(f  —  s)  by  other  finite  measures  of  past 
strain  results  in  a  different  class  of  materials  and  in 
particular  a  different  choice  for  /  in  (17). 

On  setting  £,  =  £j  in  (19)  and  using  (10),  we  see 
that 


!£(t),  if  t  >t\ 
e2,  if  0  <  t  <  f ;  (16) 

0,  if  T  <  0 ; 

where  e( x)  is  any  continuous  function  of  r  such  that 
e(t)  =  £| .  Let  U( t)  denote  the  acceleration  wave 
speed  at  time  t  >  t  for  the  strain  history  (16)  ahead 
of  the  wave.  Then  on  taking  the  limit  as  T  approaches 
t  from  above,  we  obtain  (15).  That  is,  (15)  gives 
the  acceleration  wave  speed  in  the  state  immediately 
following  the  second  strain  jump  (at  time  t)  for  the 
strain  history  (11). 


G(£i,f)  —  G(£j,0) 


gff(£l,f)^g;(£l) 

/(£  1) 


(22) 


Now  consider  an  acceleration  wave  in  a  material 
governed  by  the  integral  law  (17),  with  the  material 
ahead  of  the  wave  undergoing  the  stress  relaxation 
test  (8).  Then  from  (14),,  (19)  and  (22),  the  wave 
speed  at  time  t  >  0  is  given  by 


p0u2(t) 


di<>R{evt) 

+  y^y  [^(et)-^(et-?)]- 


(23) 
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Note  that  under  the  given  assumptions,  <y  (£j,l)  is 
the  stress  at  the  wave  front.  On  taking  the  limit  as 
t  — >  °°  and  using  (9)2,  we  obtain  the  speed  of  a 
wave  propagating  into  a  region  which  has  been  in 
equilibrium  for  all  time  at  strain  £j  and  stress  (fij ) : 

'V'H s?+7i|y <24) 

For  the  special  case  where  /  is  given  by  (20) ,  this  re¬ 
duces  to  the  formula  (1)  of  Nunziato  et  al.  [1,  §21], 
They  used  this  to  calculate  the  speed  of  expansive  ac¬ 
celeration  waves  propagating  into  deformed  regions 
in  equilibrium,  the  deformation  having  been  induced 
by  the  passage  of  a  steady  shock  wave.  The  strain 
history  resulting  from  the  passage  of  a  steady  shock 
is  only  approximately  given  by  (8),  but  due  to  fading 
memory  this  approximation  leads  to  small  errors. 


PIPKIN-ROGERS  MATERIALS 


Next  we  consider  single  integral  laws  of  the  form 

f°°  (25) 

+  h(e(t))-  J  G' (e(t  —  s),s)  ds , 


where  G'(e2is)  =  J^G(£2,s),  G(0,s)  =  0,  and 
h( 0)  =  1 .  This  is  the  special  case  of  (3)  with 


y'(el,e2,s)=h(e1)-G'{e2,s).  (26) 


Condition  (5)1  is  satisfied,  and  (5)2  holds  if  we  take 


&(£l,£2,s)  =  h(£ J  [G(£2,s)  -G(e2,0)] 

+  °)(ei) ■ 


(27) 


Pipkin  and  Rogers  [4]  considered  the  3-D  single 
integral  law 


poo 

S(f)=G(E(t),0)+y  G'(E  (t-s),s)ds,  (28) 

where  G(0,.v)  =  0,  so  that  the  first  term  on  the  right 
represents  the  instantaneous  elastic  response.  For  the 
1-D  case  considered  here,  (28)  reduces  to  (25)  with 
h(s)  =  1  —  £  =  Ap  this  term  is  a  consequence  of  con¬ 
version  from  the  second  to  the  first  Piola-Kirchhoff 
stress  measure  (see  (2)).  Replacement  of  S  in  (28) 


with  other  (Lagrangean)  stress  tensors  would  result 
in  different  functional  forms  for  h  in  (25). 

Now  consider  an  acceleration  wave  in  a  material 
governed  by  the  integral  law  (25),  with  no  restric¬ 
tions  on  the  (uniaxial)  strain  history  ahead  of  the 
wave.  From  (13)  and  (26),  we  see  that  pQU2(t)  is 
given  by  cr'(£(f))  +h'(E(t))  /0”G'(£(f  —  s),s)  ds. 
Then  on  solving  (25)  for  this  integral  and  substitut¬ 
ing  the  result  into  the  above  expression,  we  obtain 


P0U2(t)  =  a/(e(t)) 
h'(e(t )) 
h(£{t)) 


a;(£(f))-cj(f)], 


(29) 


with  £(t)  and  o(t)  the  strain  and  stress  at  the  wave 
front.  When  h(s)  =  1  —  £,  the  term  h'(£)/h(£)  re¬ 
duces  to  1/(1  —  £)  =  1  /Aj  .  For  a  stress  relaxation 
test  (8)  ahead  of  the  wave  and  t  >  0  ,  (29)  reduces  to 


,  do 


-“»(*!•')]  •  <30> 

Finally,  we  note  that  (29)  and  (30)  remain  valid  if 
(25)  is  replaced  by  the  more  general  relation 

a(t)  =  o )(e(t))  +h(e(t))-  y/(s(t-s)),  (31) 

i>0 


given  appropriate  restrictions  on  the  functional  i jf.  In 
particular,  it  is  assumed  that  t f/  does  not  depend  on 
the  current  strain  e(t).  This  includes  the  case  where 
y/  is  given  by  multiple  hereditary  integrals  of  the  type 
considered  by  Green  and  Rivlin  [6], 
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